February 14, 2013 2:18 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



NUCLEOSYNTHESIS IN THE GAMMA-RAY BURST ACCRETION 
DISKS AND ASSOCIATED OUTFLOWS 

INDRANI BANERJEE, BANIBRATA MUKHOPADHYAY 

Department of Physics, Indian Institute of Science, 
Bangalore 560012 
E-mail: indrani@physics.iisc.ernet.in , bm@physics.iisc.ernet.in 

We investigate nucleosynthesis inside the gamma-ray burst (GRB) accretion disks formed 
by the Type II collapsars and outflows launched from these disks. We deal with accretion 
disks having relatively low accretion rates: O.OOIMqs"^ Aif < O.OIA/qs"^ and hence 
they are predominantly advection dominated. We report the synthesis of several unusual 
nuclei like "^^P, ^^K, *^Sc, ^^Cl and various isotopes of titanium, vanadium, chromium, 
manganese and copper in the disk. We also conflrm that isotopes of iron, cobalt, nickel, 
argon, calcium, sulphur and silicon get synthesized in the disk, as shown by previous 
authors. Much of these heavy elements thus synthesized are ejected from the disk and 
survive in the outflows. Indeed, emission lines of many of these heavy elements have been 
observed in the X-ray afterglows of several GRBs. 
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1. Introduction 

The coUapsar model is the most promising theoretical model explaining the long- 
duration gamma-ray bursts (GRBs) and the supernova (SN)— gamma-ray burst con- 
nections.^ We plan to investigate the nucleosynthesis in the accretion disks formed 
by the Type II collapsars^ where the accretion rate (M) is: Q.QQ\MqS~^ ^ ^ < 
0.01Mqs~^ , when Mq indicates solar mass, as this regime is the site for the syn- 
thesis of heavy elements. Neutrino cooling becomes important in the inner disk 
where the temperature and density are higher. We report the synthesis of various 
elements like '^^P, "^^K, ^^Sc, '^^Cl and several uncommon isotopes of titanium, vana- 
dium, chromium, manganese and copper in the disk for the first time, to the best of 
our knowledge, apart from isotopes of iron, cobalt, nickel, calcium, argon, sulphur 
and silicon which have already been reported earlier.t2l 

We also consider nucleosynthesis in the outflows from these disks using an adia- 
batic one-dimensional outflow modeP' and report that many of the heavy elements 
thus synthesized in the disk do survive in the outflow. Additionally, depending on 
whether significant '^^Ni is synthesized in the outflow, it can be predicted whether 
the outflow will lead to a supernova explosion or not. 



2. Disk and Outflow Models 

We adopt height-averaged equations based on a pseudo-Newtonian framework.Sl 
The accretion disk formed in a Type II collapsar is modelled within the framework 
suggested by previous authors]^ where the electron degeneracy pressure and the 
evolving neutron to proton ratio are appropriately calculated. 
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To investigate nucleosynthesis in the outflow, we adopt a simple, one- 
dimensional, spherically expanding and adiabatic outflow modelP 

3. Nucleosynthesis inside accretion disks and outflows 

We use He-rich and Si-rich abundances as the initial conditions of nucleosynthesis at 
the outer disk. Well tested nuclear network codej^has been used. We have modified 
this code further by increasing the nuclear network and including reaction rates 
from the JINA Reaclib Database, https://groups.nscl. msu.edu/jina/reaclib/db/.'^ 

We also consider outflow from various radii of ejection, R^j, with Rej < 200Rg, 
and evaluate the abundance evolution in the outflow assuming the initial composi- 
tion the same as in the accretion disk at R^j . 

Figure 1 (a) models the abundance evolution around a SMq Schwarzschild black 
hole accreting at AI — 0.001Mqs~^ , with the viscosity parameter a = 0.01 and the 
composition of the accreting gas at the outer disk similar to the pre-SN He-rich layer. 
It depicts that the disk comprises of several zones characterized by the dominant 
elements. The region ^ 1000 — 300i?g, Rg being Schwarzschild radius, is mainly the 
^''Ca, ^''Ti and "^^Cr rich zone. This is because unburnt "^^Ar undergoes a-capture 
reaction to give rise to '*°Ca through ■^^Ar(Q!, 7)^°Ca which further undergoes a- 
capture reactions to give rise to *^Ti and "^^Cr. 

Inside this region, the temperature and density in the disk favor complete photo- 
disintegration of '^''^Ti and '*^Cr releasing a-particles and resulting in the formation 
of ''°Ca, ^'^Ar, ^^S and ^^Si, as is evident from Fig. 1(a). Subsequently, ^^Si and 
^^S start burning, which favors formation of iron-group elements via photodisinte- 
gration rearrangement reactions.'^ Therefore, in the range ^ 300 — 8QRg, there is 
a zone overabundant in ^^Ni, ^''^Fe, '^^S and ^^Si. Inside this zone, all the heavy 
elements photodisintegrate to '^He, neutron and proton. 

On taking the same disk and black hole parameters but considering Si-rich abun- 
dance at the outer disk we note that the disk has a huge zone rich in ^*Si and "^^S 
extending upto radius ^ 250i?g. Inside this radius, silicon burning commences and 
soon the disk becomes rich in '^'*Fe, ^^Ni and ^^Ni. Inside 7QRg, all the heavy ele- 
ments again get photodisintegrated to ^He, protons and neutrons. Another remark- 
able feature in the He-rich and Si-rich disks is that inside ^ 100i?g, the abundances 
of various elements start becoming almost identical as if once threshold density and 
temperature are achieved, the nuclear reactions forget their initial abundances and 
follow only the underlying disk hydrodynamics. 

On increasing M ten times we find that the individual zones shift outward 
retaining similar composition as is in the low M cases described above. 

Next we consider the outflow from 40i?g of the He-rich disk with M = 
O.OOIMqs"^, as described above. We show that ^^Ni is abundantly synthesized in 
the outflow along with isotopes of copper and zinc, as shown in Fig. 1(b). Synthesis 
of ^^Ni signifies that the outflow will lead to a supernova explosion. Figure 1(b) 
also depicts that on changing the initial velocity of ejection the final abundances 
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Fig. 1. (a) Zones characterized by dominant elements in the disk with He- rich abundance at the 
outer disk, (b) Abundance evolution in the outflow from R^j ~ AORg of the disk in (a) , solid lines 
correspond to the higher velocity of ejection and other lines to lower velocity, in each set. 



of the nucleosynthesis products change significantly. More ^^Ni is synthesized when 
the velocity of ejection is low (see Fig. 1(b)) because then the temperature drops 
slowly in the ejecta which facilitates greater recombination of alphas to nickel. 



4. Conclusions 

We report for the first time, to the best of our knowledge, that several unusual nu- 
clei like ^^P, ^®K, ^^Sc, ^^Cl and various uncommon isotopes of titanium, vanadium, 
chromium, manganese and copper get synthesized in the disk apart from isotopes 
of iron, cobalt, nickel, silicon, sulphur, argon and calcium. Several zones, character- 
ized by dominant elements, are formed in the disk. The elemental distribution in the 
disk starts appearing identical once threshold density and temperature are reached. 
The respective zones shift outwards on increasing the mass accretion rate. Several 
of these heavy elements survive in the outflow from these disks, and when ^^Ni is 
abundantly synthesized in the outflow, there is always a supernova explosion. 
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